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Abstract 

In supersymmetric models where the magnitude of the GUT scale gaugino mass parameter 
M3 is suppressed relative to M\ and M2, the lightest neutralino can be a mixed higgsino- 
bino state with a thermal relic abundance in agreement with the WMAP central value 
for OcDM^ 2 and consistent with all other phenomenological constraints. In these models, 
the gluino can be as light as 200 GeV without conflicting with the LEP2 bounds on the 
chargino mass. Thus, gluino pair production can be accessible at the Fermilab Tevatron 
at high rates. In this framework, gluinos decay radiatively with a large branching fraction 
to a gluon plus a neutralino. We find that experiments at the Fermilab Tevatron, with 5 
fb _1 of integrated luminosity, will be sensitive to gg production in the mg ~ 200 — 350 
GeV range via the multi-jet +i?™ lss and multi-jet +£ + £~ + i?™ lss channels at the 5cr level, 
while trilepton signatures are expected to be below this level of detectability. Dilepton 
mass edges from both Z2 and Z3 decays may be measurable in the dilepton + multi-jet 
+Ef iss channel. 
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1 Introduction and motivation 



The determination of the average density of cosmological cold dark matter (CDM) |]Q 

n CBM h 2 = o.minis (2<r), (i) 

imposes a stringent constraint on any beyond the Standard Model framework featuring a weakly 
interacting massive particle stable on cosmological time-scales. 1 In particular, ((TJ) poses a severe 
constraint on i?-parity conserving supersymmetry (SUSY) models where the lightest neutralino 
(Zi) is the lightest supersymmetric particle (LSP) [2]. 

Although it is possible to reconcile the value of f2cDM^ 2 determined by the WMAP team £Q 
with the thermal neutralino relic abundance fi^/i 2 expected in the framework of the minimal 
supergravity (mSUGRA) model [3], agreement with is obtained only within narrow regions, 
most of which are close to the boundary of the allowed parameter space. While the smallness of 
these regions reflects the impressive precision achieved in the determination of Ocdm, the fact 
that they lie close to phenomenologically constrained portions of the parameter space reflects a 
general result in the mSUGRA setup: except in the case where sparticles are light (the so-called 
bulk region), f2cDM^ 2 is considerably smaller than the typical mSUGRA expectation for VL-^Ji 2 . 
Special neutralino annihilation mechanisms can, however, be operative in the Early Universe, 
enhancing the LSP pair annihilation rate and consequently suppressing its relic abundance 
to acceptable values. In mSUGRA, instances of such mechanisms are resonant neutralino 
annihilations through s-channel Higgs exchange diagrams j3], the edges of parameter space 
where the LSP co-annihilates [S] with either a light stau [Bj or a light stop [7], or where |//| 
is small enough so that the LSP features a substantial higgsino component (the hyperbolic 
branch/focus point (HB/FP) region) jS]. Several groups have examined the signals expected 
in collider experiments, as well as via direct and indirect searches for neutralino dark matter 
in underground detectors, assuming that the parameters are in one of these WMAP-allowed 
regions of the mSUGRA model. 

Motivated by the fact that the correlations between the WMAP measurement and expec- 
tations in other experiments may be mo del- dependent, there have been a number of recent 
studies that have relaxed the universality assumption, that is the hallmark of the mSUGRA 
framework. Allowing for non-universal Higgs boson mass (NUHM) parameters allows for an 
extended region of parameter space where resonant annihilations occur or/and where \y\ is 
sufficiently small jH], while non-universality in the SU(2) and U(l) GUT scale gaugino mass 
parameters allows agreement between and ^cdm either via an enhanced wino fraction in 
the LSP PU]; or via bino-wino coannihilation [TJ. These extended scenarios can be distin- 
guished from one another, as well as from the minimal mSUGRA scenario, because they give 
rise to different outcomes for collider signals and for the anticipated detection rates at dark 
matter search experiments. 

Another option to obtain a consistent thermal neutralino relic density is to reduce the 
magnitude of the GUT scale SU (3) gaugino mass M 3 relative to the magnitude of the SU (2) 
and U(l) gaugino masses ^21 El (the so-called low |M 3 | dark matter model (LM3DM)). As 

1 We quote the value obtained by the WMAP collaboration by combining their data with that from the Sloan 
Digital Sky Survey. 
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explained in Ref. a lowered relic density occurs because a smaller value for |M 3 | also 
induces lower values for the squark masses and the trilinear scalar couplings via the running 
dictated by the (coupled) renormalization group equations (RGEs). The RGE running also 
yields a suppression in the absolute size of the soft breaking Higgs mass squared parameter 
\m 2 H |, which, in turn, lowers the magnitude of the weak scale superpotential mass parameter 

(fixed by the value of Mz), so that the lightest neutralino develops a significant higgsino 
component, giving rise to mixed higgsino dark matter (MHDM) 2 . Agreement with WMAP 
is then obtained because the neutralino annihilation rate is enhanced by larger annihilation 
amplitudes into gauge and Higgs boson pairs, and co-annihilation with the lightest chargino 
and the next-to-lightest neutralino further suppress the final LSP relic density. 

In the LM3DM scenario, we generically expect the ratio of the gluino to lightest chargino 
mass nig : to be smaller than the corresponding value ~ 3 — 3.5 expected in models with 
universal GUT scale gaugino masses and large This ratio is important when comparing 
collider searches for sparticles with LEP and the Fermilab Tevatron. Assuming that ~ m z 1 
is not too small and that rrty > 200 GeV, consistency with LEP2 experiments requires > 
103.5 GeV. For models with gaugino mass unification and large \fi\, this bound implies that 
gluinos must have mass greater than ~ 300 — 350 GeV. Such large mass gluinos are difficult 
to search for at the Tevatron, as their production cross section is rapidly suppressed with 
increasing masses. 

In the LM3DM model, instead, relatively light gluinos (values of the gluino mass as low 
as nig ~ 200 GeV would be consistent with the LEP2 constraints) can be copiously produced 
in hadronic collisions, and the currently operating Fermilab Tevatron is the obvious facility 
to search for these new matter states. To-date, experiments at the Tevatron have searched 
for gluinos in their multi-jet plus E™ 1SS data sample, and exclude gluinos lighter than roughly 
200 GeV, irrespective of the squark mass, from their analysis of the Run 1 data [T^]. Very 
recently, the D0 collaboration, from an analysis of 310 fb _1 of data has obtained a new limit JH] 
of m g > 233 GeV. Unlike the multi-lepton plus jets plus E™ 1SS analyses based on cascade decays 
of gluinos, inclusive E™ ss analyses are largely independent of the details of the spectrum in the 
electroweak "-mo" sector 3 . 

Within any framework with unification of gaugino masses, m g ~ (3 — 3.5)771^, and the 

published limits from CDF and D0 are pre-empted by the LEP limit ~ 103 GeV on the 
chargino mass. Within the LM3DM scenario, instead, the gluino is relatively light, and the 
impact of the LEP chargino limit on the Tevatron gluino search is clearly reduced, so that 
it is possible that data from Tevatron experiments may probe a range of the LM3DM model 
parameter space not accessible to LEP2, either in the current data sample, or in the data 
sample expected to be accumulated at the Tevatron before the LHC completes about a year of 
operation. 

In this study we explore the prospects for detection of gluino pair production within the 
framework of the LM3DM scenario. The remainder of the paper is organized as follows. In Sec. 

2 Although the gluino mass is reduced from its usual value, the gluino - LSP mass splitting is still large so 
that gluino co-annihilation can be safely neglected in the evaluation of the Z\ relic density ^1] . 

3 These analyses are not completely independent of chargino and heavier neutralino masses because the 
transverse momenta of the Z\ LSPs, and hence the E^ lss spectrum, does depend on the cascade decay patterns. 
Moreover, sometimes a lepton veto is also imposed on the SUSY signal. 
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2 we discuss the parameter space and sparticle mass spectra expected in the LM3DM model. 
In Sec. 3, we discuss signal rates and backgrounds for gluino pair discovery in the jets +E™ 1SS 
channel. In Sec. 4, we show that discovery in the clean trilepton -\-E™ lss channel is unlikely. 
In Sec. 5, we show that detection in the dilepton plus multi-jet +EJp lss is a viable possibility, 
and that the associated m(£ + £~) distribution can give the characteristic mass edges indicative 
of the — m^, and, possibly, also of the — mass difference. Finally, we summarize 
our results in Sec. 6. 



2 The low IM3I dark matter model 

The low IM3I dark matter model differs from mSUGRA only in that the GUT scale gluino mass 
parameter M3 needs not be equal to m^j = M\ = Mi- The parameter space of this model is 
thus given by, 

m ,m 1/2 ,M 3 , A , tan/3, sign(fi), (2) 

where mu2 is taken to be positive without loss of generality, but M 3 can take either sign. 
For any set of values for the parameters in (J2|l . we can vary r% = M3/7TI1/2 so as to increase 
the higgsino content of the LSP and to drive the LSP annihilation rate to yield a relic LSP 
density fi^ h 2 in agreement with (0). In order to get small enough, we must "slow down" 
the RG evolution of rn 2 H from its GUT scale value of m\ to a negative value at the weak scale 
- remember that m 2 Hu (weak) ~ — fi 2 as long as tan (3 is not very small - which, in turn, requires 
a smaller value of X t = + mf + m 2 Hu + A 2 than in mSUGRA. Since gauge coupling effects 
always increase squark mass parameters as they evolve from the GUT scale down to the weak 
scale, and since the large SU(3) gauge coupling contributes dominantly to this increase, smaller 
values of X t are obtained by choosing |M 3 (GUT)| to be smaller than its mSUGRA value of 
mi/2- 

We provide a panorama of the LM3DM scenario in Fig. Q where we show contours of fixed 
gluino mass in the tuq — m 1( / 2 plane with tan/3 = 10, Aq = and fi > 0, and where at each 
point in this plane r 3 has been chosen to obtain the central value given in (JIJ for the LSP 
relic density. We use Isajet v7.74 for sparticle mass calculations ^7]. The grey (red) region 
is excluded because either electroweak symmetry in not properly broken or the LSP becomes 
charged or colored. The black (blue) region is excluded by the LEP2 negative search results 
for charginos. The wiggles in the plot curves reflect numerical issues related to the precision 
with which we require the neutralino relic abundance to saturate the WMAP central value for 
the CDM abundance (JTJ, and also any numerical instabilities in the code for the determination 
of /1 as a function of r 3 . On the extreme left of the plot where the gluino mass contours dive, 
the Z\ is dominantly a bino since (due to light sleptons) the r% value there need not deviate 
severely from ~ 1. As we move to larger values of mo at fixed m^, much smaller values of r% 
are needed for the neutralino relic abundance to match the CDM density in (0), and we step 
into the MHDM region which, as explained above, also features a small value of m g . Indeed 
we see that for mo ~ 1 TeV, the gluino could be lighter than even 200 GeV in a region of 
parameter space unconstrained by the negative results of sparticle searches at LEP2. 

As an example of the relation between sparticle masses in this region of parameter space, we 
show in Fig. |2]the value of m,g, together with the chargino and neutralino masses (the sfermions 
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LM3DM: -M 3 < m V2 , tan$=10, A =0, n >0, m, =175 GeV 
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Figure 1: Contours of nig in the tuq vs. m\ji plane for the LM3DM model, where M^{Mgut) has 
been set, at every point of the parameter space, to the value such that Vt~ h 2 = 0.11. We take Aq = 0, 
tan/3 = 10, mt = 175 GeV, and consider M 3 < in the left panel and M 3 > in the right panel. The 
grey (red) regions are excluded because either the electroweak symmetry cannot be correctly broken, or 
because the LSP is charged. The black (blue) shaded regions are excluded by the LEP2 bound on the 
chargino mass. 



are too heavy to be accessible at the Tevatron) versus m 1( / 2 for the slice of the plane in Fig.^at 
fixed m = 1500 GeV. This ijiq value is representative of the range needed for which |M 3 (GUT) | 
has to be significantly reduced from its mSUGRA value in order to obtain agreement with the 
observed value of ^cdm^- 2 - While in mSUGRA one expects the masses rrtg : : to be in 
the ratio ~ 7 : 2 : 1, we find here that with MHDM, the typical ratio is rather ~ 2.5 : 1.5 : 1, so 
that not only is the m g — mass gap reduced, but the — mass gap is suppressed as 
well. Another noteworthy feature is that because of the smallness of there is sizable mixing 
between gauginos and higgsinos resulting in three relatively light neutralinos, while the heavy 
chargino and the heaviest neutralino (which are dominantly wino-like) are considerably split 
from their lighter siblings. While all the masses increase steadily with mi/2, for the M3 > 
curves (solid lines) we see sharp glitches at very low m\/2 ~ 270 GeV where < Mw' for 

mi/2 < 270 GeV, very low values of r 3 are needed since Z^Zi — > W + W~ annihilation in the 
early universe becomes kinematically suppressed. There are similar glitches for negative M 3 
(dashed lines), but these occur for mu2 values excluded by the LEP2 constraints, and are not 
seen in the figure because we terminate the curves on the left when the chargino mass falls 
below its LEP2 limit. 

Prior to discussing whether Tevatron experiments can probe supersymmetry in this region 
of parameter space we need to study the decay patterns of the gluino and of its daughter 
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LM3DM: m a =1500GeV, tan$=10, A =0, \i >0, m t =175GeV 
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Figure 2: Gluino, chargino and neutralino masses versus mi/2 for the LM3DM model where 
M%(Mgut) has been lowered at every point to obtain ^z_h 2 = 0.11. We take here Aq = 0, tan/3 = 10, 
mt = 175 GeV and mo = 1500 GeV. The solid curves correspond to M3 > 0, while the dashed curves 
to M3 < 0. We cut the curves on the left when the chargino mass falls below its LEP2 bound. 

sparticles. For reasons detailed in Ref. ^3], the radiative decays g — > gZi dominate for the 
gluino masses of interest at the Tevatron. In the upper frames of Fig. El we show the branching 
ratio for these various radiative decays of the gluino for M 3 > (left frame) and M 3 < (right 
frame), together with that for the sum of all its three-body decays (labeled 3), versus m\/2- We 
adopt here the same parameter set as in Fig. |21 As in the preceding figures, we set M 3 so that 
Sl^h 2 = 0.11, the WMAP central value for the CDM abundance pQ. We see that - depending 
on the sign of M 3 - gluinos lighter than ~ 420 — 475 GeV dominantly decay radiatively. For 
small to medium values of tan /3, where bottom quark Yukawa couplings can be neglected, the 
partial width for the various radiative decays is mainly governed by the H u content of the 
neutralino f3J EI], and accounts for the ordering of the branching fractions for these decays. 
The sharp rise in the branching fraction for the three body decays is due to the opening up of 
decays to the wino-like Z 4 and W2, both of which have large SU(2) gauge couplings to q^. when 
these modes are not phase space suppressed, they rapidly dominate the decay width. Note that 
unlike the tree-level decay, the radiative decay to the dominantly wino-like Z 4 is dynamically 
suppressed because the higgsino component of the wino-like state is always small. Although we 
have shown the results for the particular choice of mo = 1.5 TeV, we have checked that these 
results are qualitatively unaltered for m values in the range between 1 and 2 TeV. 

For large values of tan j3, the value of r 3 required to get fl^ h 2 = 0.11 is larg Gr^ clS £L result, the 

gluino is relatively heavier than Z4 or W2, and its tree- level decays become dominant. However, 
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LM3DM: m Q =1500GeV, tan§=10, A Q =0, \i >0, m { =175GeV 
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Figure 3: The curves (labeled Zi) in the upper panels show the branching fractions for the radiative 
decays g — > Zig of the gluino versus m^j for LM3DM model where M^{Mqut) has been adjusted at 
every point to attain mixed higgsino dark matter with Q^h 2 = 0.11, while the curve labeled 3 denotes 
the corresponding branching fraction for the sum of all three body decays of the gluino. The (barely 
visible) unlabeled lowest curves in the upper panels indicate B(g — > Z±g) which is find to lie always 
below the percent level. The lowerjrames show the total leptonic branching fractions for the decays, 
Z2 — ► ££ + X, Z3 — > ti + X and W\ — > £ + X, adding in all possible decay chains for the particular 
lepton topology. The left (right) frames are for M3 > (M3 <0). Everywhere, we fix mQ = 1500 GeV, 
Aq = 0, tan/? = 10 and mt = 175. 

bearing in mind the possibility that the cosmological dark matter can very well consist of several 
components, lower values of r 3 (which would lead to fig- < ficDM) are not excluded. Within 
the LM3DM framework, Tevatron experiments can, and should, search for gluinos in this large 
tan/3 portion of the parameter space since it has not been excluded by LEP2 searches. By 
the same token, if dark matter has several components, it is possible that the gluinos are even 
lighter than what we obtain here (see, e.g. Fig. |2J), and the resulting LEP2 excluded region 
could as well be smaller than what we show. 

In the lower frames of Fig. El we show the cumulative leptonic branching fractions for the 
daughter neutralinos and charginos versus mm. For W\ and Z 2 , this is simply the usual 
branching ratio B{W\ — > tvZ\) and B{Zi — > ££Zi), but for Z3 the two leptons can come from 
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either its primary decay, or from the leptonic decays of daughter neutralinos 4 . The branching 
fractions shown in these lower frames can be used in conjunction with those in the upper frames 
and the gluino production cross sections to estimate cross sections (before any cuts) for various 
multi-lepton topologies in di-jet events at the Tevatron. 

3 Search in the jets +E™ SS channel 

In this section, we examine whether the Fermilab Tevatron can detect gluino pair production 
in the LM3DM model in the multi-jet -\-E™ lss mode, assuming 5 fb _1 of integrated luminosity 
that is projected to be accumulated by each experiment at the Tevatron. We generate signal 
and background events using Isajet 7.74, with a toy detector simulation containing hadronic 
calorimetry ranging out to \rj\ < 4, with cell size Arj x A0 = 0.1 x 0.262. We adopt hadronic 
smearing of AE = 0.7 j\[E and EM smearing of AE = 0.15/v 7 ^. We adopt the Isajet GET JET 
jet finding algorithm, requiring jets in a cone size of AR = 0.5 with > 15 GeV. Jets are 
ordered from highest E T (ji) to lowest E T . Leptons within \r)f\ < 2.5 {I = e, /i) are classified 
as isolated if Pt{£) > 5 GeV and a cone of AR = 0.4 about the lepton direction contains 
Et < 2 GeV. Finally, if a jet with \rjj\ < 2 has a 5-hadron with > 15 GeV within 
AR < 0.5, it is tagged as a 6-jet with an efficiency of 50%. 

To find optimal cuts, we generated 100K signal events for the case where mi/2 = 300 
GeV, m = 1500 GeV, A = 0, tan (3 = 10 and \i > 0. For this point, M 3 = 79.69 GeV 
yields ^^_h 2 = 0.12. We have also generated SM background event samples from W + jets 
production, Z + jets production, ti production and vector boson pair production 5 . The W or 
Z + jets sample uses QCD matrix elements for the primary parton emission, while subsequent 
emissions are generated from the parton shower. We adopt a set of cuts similar to those used by 
the D0 collaboration in Ref. [T^j, but optimize the E™ 1SS and Ht cut values for this framework. 
Our final set of cuts are listed in Table ^ where we divide the signal topologies into > 2- 
jets+.E'™ lss , > 3-jets+i?™ lss and > 4-jets+i?™ lss , while vetoing isolated leptons. The constituent 
background rates from the major background sources are listed in Tabled From these rates, we 
can compute the signal observability level needed for a given integrated luminosity, using the 
following criteria: (1) the statistical significance S/\/B > 5cr, (2) S/B > 25%, and (3) S > 10 
events. 

Our results for the SUSY reach of the Tevatron within the LM3DM framework are shown 
in Fig. Inversus mm for the same parameter choices as in Fig. |21 Assuming an integrated 
luminosity of 5 fb _1 , we have checked that the reach in each of the three n-jet + E™ 1SS event 
topologies is limited by the 5a criterion. The minimum cross section for observability of the 
signal is shown by the dashed horizontal line, while the signal is indicated by the solid (dashed) 
curve for M 3 > (M 3 < 0) for (a) > 2-jets+££ liss events, (6) > 3-jets+£™ iss events and 
(c) > 4-jets+E'™ lss events. We see in each of frames (a)-(c) that the 5 fb _1 reach extends out 
to mi/2 ~ 330 — 340 GeV, corresponding to a reach in rrig according to Fig. [21 of ~ 320 GeV. 
Within the LM3DM framework, this corresponds to a reach in ~ 170 GeV, and thus 

4 In principle, there could be contributions from Z 3 — > W^W^ 1 , but these decays are kinematically inaccessible 
over the entire parameter space range shown in the plot. 

5 We do not estimate QCD backgrounds which, we assume, are negligible after the cuts described below |16| . 
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cut 


2j + 


3j + £^ 1SS 


Aj + E^ 1SS 


A0(j 1 ,j 2 )<165° 


yes 


yes 


yes 


isol. lep. veto 


yes 


yes 


yes 


rij 


> 2 


> 3 


> 4 


fe,l<0.8 


3i, 32 


31J2J3 


jl, 32, j 3, 34 


80° < A<^>(£^ iss ,ji) < 150° 


yes 


yes 


yes 




50° - 150° 


50° - 150° 


60° - 150° 


771 miss 


> 120 GeV 


> 100 GeV 


> 75 GeV 


Ht 


> 220 GeV 


> 150 GeV 





Table 1: Cuts used for the analysis of multi-jet +£ , ™ 1SS signatures in the LM3DM model. 



BG 


2j + E^ ms 


3j + E!£ 1SS 


4j + E^ ms 


tt(175) 


6.6 ±0.3 


12.3 ±0.5 


14.9 ±0.6 


W + jets 


8.9 ± 1.4 


15.5 ± 1.9 


12.1 ± 1.7 


Z + jets 


11.0 ±0.7 


17.2 ±0.9 


9.0 ±0.7 


total 


26.5 


45.1 


36.0 



Table 2: STW backgrounds in fb a/ter ra£s /istec? in TableUlfor the multi-jet +£ , ™ 1SS signatures in the 
LM3DM model. 

extends well beyond that of LEP2 experiments. 

4 Search in the trilepton +EJf lss channel 

We have also examined the reach of the Fermilab Tevatron in the much touted inclusive trilepton 
channel ^H] where the leptons arise from the decays of charginos and neutralinos produced via 
pp — > W\Zi + X, or via cascade decays of gluinos. Since, as discussed above, the mass gap 
between Wi/Z 2 and the Zi LSP is not large, we expect the lepton spectra to be relatively 
soft. Hence, for this study, we adopt the soft lepton cuts SC2 introduced in the first paper of 
Ref. [20], where the background was found to be 1.05 fb. The reach in the inclusive trilepton 
channel is shown in Fig. Efd) where we see that signal is always below the 5<r observability 
level. This is, in part, due to the fact that the kinematically favored W\Z 2 ^ production now 
dominantly occurs via the weak isodoublet higgsino components of the chargino and neutralino 
which have a smaller coupling (than the weak iso-triplet coupling characteristic of the mSUGRA 
framework) to the Z boson. We conclude that in the case of the LM3DM model, the best search 
channel is the multi-jets +EJp lss channels. 

5 Search in the jets + OS-dilepton +E™ SS channel 

The relatively low value of is the characteristic feature of the LM3DM model. As a result, 
three (rather than two) neutralinos tend to be relatively light and mixed, whereas gaugino- 
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Figure 4: The SUSY reach of the Fermilab Tevatron with 5 fb^ 1 of integrated luminosity in the 
multi-jet +£vp lss channel for (a) di-jet events, (b) tri-jet events, (c four jet events, and (d) trilepton 
events for LM3DM model where \M^(Mgut)\ has been adjusted at every point to get VL~ h 2 = 0.11. 
We fix A() = 0, tan (3 = 10, nit = 175 GeV and mo = 1500 GeV, a slice representative of the LM3DM 
parameter space under investigation here. The solid curve is for M3 > 0, while the dashed curve 
corresponds to M3 < 0. 

higgsino mixing increases the masses of the heavier chargino and the heaviest neutralino. It is, 
therefore, reasonable to ask whether it is possible to identify their production via the cascade 
decays of gluinos at the Tevatron. We are thus led to investigate the observability of the signal 
in the multi-jet + opposite sign (OS) dilepton + E™ lss channel, where the leptons have the same 
flavor. This channel is of special importance since it has been long known that the dilepton 
invariant mass spectrum from Zi — > Zi+££ contains a kinematic cut-off at m-x —m^ . The mass 

Zji Zi\ 

edge(s), if visible, can serve as the starting point for reconstructing sparticle cascade decays, 
and for obtaining information on sparticle masses [2~T] . 

Toward this end, we examine the signal in the multi-jet +££ + E™ 1SS channel, where £ = e 
or /i. We extract signal events containing two opposite-sign/same flavor isolated leptons plus 
jets plus missing transverse energy, and compare the signal with SM backgrounds from ti 
production, Z — > rf + jets production and vector boson pair production Z°Z° and 

W ± Z° production). By requiring hard missing (E™ ss > 75 GeV), we reject much of the 
background from Z° production, while by requiring a veto of events with a tagged fe-jet we reject 
much of the ti background with hardly any loss of signal. Finally, requiring at least 2 jets in 
the events improves the statistical significance of the signal. The surviving background rates in 
fb, along with signal in the LM3DM framework for = 300 GeV and other parameters as in 
Fig. 12} are listed in Table El The corresponding reach in the > 2-]ets+££ + E™ 1SS channel is once 
again governed by the 5a criterion, and is shown in Fig. Inversus mi/2, with other parameters 
as in Fig. |21 for 5 fb" 1 of integrated luminosity. We see that it extends out to m 1( / 2 ~ 310 — 320 
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BG 



2j 



/rimiss 



tt(175) 11.6 ±0.5 

Z —> rf + jets 5.6 ± 0.5 

WW, WZ, ZZ 7.6 ±0.6 

total 24.8 

signal m 1/2 = 300 GeV 21.4 ± 0.6 



Table 3: SM backgrounds and signal for mj m = 300 GeV in fb after cuts listed in text for the multi-jet 
+££ + Ef iss signatures in the LM3DM model. 



LM3DM: m g =1500GeV, tan$=10, A Q =0, \i >0, m f =175GeV 

,2 




250 275 300 325 350 375 400 425 450 475 500 

m 1/2 (GeV) 



Figure 5: The SUSY reach of the Fermilab Tevatron in the > 2 jets+££ + Ej^ lss channel for the case 
of mo = 1500 GeV, Aq = 0, tan/3 = 10 and mt = 175 GeV assuming an integrated luminosity of 
5 fb' 1 . We dial M 3 (M GUT ) for each m 1/2 so that = 0.11. 



GeV, i.e. slightly lower than the reach in the multi-jet±i?™ lss channels. 

In order to examine the detectability of any dilepton mass edges, we show the opposite- 
sign / same-flavor dilepton invariant mass spectrum from the signal and background in Fig. El for 
the case of m\n = 300 GeV, and other parameters as in Fig. El The hatched distribution comes 
from the various background sources listed in Table El which includes a peak at m(££) = Mz 
from Z-pair production. The signal plus background is shown by the open histogram. In this 
case, a distinct mass edge can be seen at —m^ ~ 59 GeV. Remarkably, the mass edge from 

Z 3 — > Z\££ is also seen at at — ~ 86 GeV. This higher mass edge will be somewhat 



obscured by Z-width effects, which are not included in our simulation of ZZ production. The 
point, however, is that the value of m^ g —m^ in our study is only fortuitously close to Mz, and 
in general, it may be possible to see even the second mass edge at the Tevatron! Observation 
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10 20 30 40 50 60 70 80 90 100 

m(tr) (GeV) 

Figure 6: The spectrum of opposite sign/same flavor dilepton invariant mass in background (hatched) 
and signal-plus-background (open histogram), for the case of = 1500 GeV, my 2 = 300 GeV, 
M3 = 79.69 GeV. We also take Aq = 0, tan/3 = 10 and m t = 175 GeV. The arrows denote the 
theoretically expected positions of the corresponding mass edges. 

of this second mass edge would provide a strong hint for a small value of 

6 Summary and concluding remarks 

Within the mSUGRA model, or any other supersymmetric setup with unification of the soft 
SUSY breaking gaugino mass parameters at the GUT scale, the lower limit > 103 GeV 
from LEP2 experiments excludes gluinos with masses smaller than about 300-350 GeV, leaving 
little room for gluino searches at the Fermilab Tevatron. This is, however, a model- dependent 
conclusion, and, as already stressed elsewhere [22], Tevatron experiments should search for 
gluinos independently of the constraints from LEP2. 

We provide here a specific example, the so-called low |M 3 | dark matter framework (LM3DM), 
where the universality of the GUT scale gluino mass parameter with the corresponding SU (2) 
and hypercharge gaugino mass parameters is relaxed, while the universality of all other soft 
SUSY breaking parameters is retained, as in the mSUGRA setup. Adjusting the magnitude 
of Ms(GUT) (which can have either sign) to low values leads to SUSY spectra with relatively 
suppressed values of \fj,\, entailing, in turn, a larger LSP higgsino fraction, which can then lead 
to an LSP relic density in agreement with the observationally preferred central value of 
JIcdm^ 2 for any value of the other soft SUSY breaking parameters. 

The LM3DM framework leads to characteristic differences in the sparticle spectra from 
the usually studied frameworks with unified gaugino masses, or with anomaly-mediated SUSY 
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breaking. In particular, low values of |M 3 | imply that the rrig/m^ ratio is significantly smaller 
in the LM3DM model compared to the mSUGRA case, so that experiments at the Tevatron 
will be able to explore regions of parameter space not already ruled out by LEP2. 

The main result of the present study is the reach of the Fermilab Tevatron experiments 
within the LM3DM framework, shown in Fig. 0] and Fig The best reach is obtained in the 
inclusive multi-jet + E™ 1SS channels, while the reach in the multi-jet plus opposite sign dilepton 
channel is only slightly less effective. Assuming an integrated luminosity of 5 fb _1 , expected to 
be delivered to each experiment within the next two years of operations at the Tevatron, the 
reach extends up to mi/2 = 350 GeV which, for M3 > corresponds to m g ~ 325 GeV and 
Triyy ~ 170 GeV, significantly beyond the reach of LEP2. Combining the two experiments will 
yield an even higher reach. 

The concomitant smallness of within this framework implies that both Z2 and Z3 may be 
accessible via gluino decays, offering another interesting opportunity to Tevatron experiments, 
as illustrated in Fig. |UJ the invariant dilepton mass spectrum for events with > 2 jets + OS 
dileptons + E™ 1SS , with a veto on 6-tagged jets (to reduce the background from ti production), 
may yield mass edges from both Z 2 — > UZ\ and Z 3 — > UZ\ decays. Observation of two mass 
edges would strongly suggest a small value of 

In summary, if SUSY is realized as in the LM3DM model, a framework consistent with all 
constraints from particle physics and cosmology, experiments at the Tevatron will be able to 
probe regions of parameter space not accessible at LEP 2 before the LHC experiments turn on 
and collect data for physics analysis. We urge our colleagues on the CDF and D0 experiments 
to search for gluinos irrespective of constraints from chargino searches since these are based on 
the untested assumption of gaugino mass unification. 
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